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ABSTRACT 

We have examined methanol emission from Orion KL with of the HerschelfHlFl instrument, and detected two methanol bands centered at 524 
GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band la) is dominated by the isolated AJ=0, K=-4— »-3, v,=0 Q branch, and 
includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the AJ=0, K=7— ^6, v,=0 
Q branch transitions which are mostly blended. We have used the isolated E-type v,=0 methanol transitions to explore the physical conditions in 
the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight 
toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in 
these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure 
in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power 
of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars. 
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' 1. Introduction 

Methanol (CH3OH) is an abundant molecule in a wid e range 
of int erstellar conditions (e.g., Friberg et al. 1988; Mente n et al.l 
I Il988l) . An asymmetric-top rotor, methanol has numerous rota- 
r tional transitions from far-infrared to millimeter wavelengths. 
This marks methanol as a common "weed" molecule, and the 
number and strength of its transitions means that methanol con- 
tamination has to be minimized when analyzing transitions from 
I other molecules. However, rather than being treated solely as a 
, contaminant, methanol is also a useful molecular tracer of dense 
gas (Menten et al. 1988). Leurini et al. (2004) further explore the 
properties of methanol emission in dense molecular clouds and 
■ illustrate how methanol transitions at sub-millimeter and mil- 
' limeter wavelengths are sensitive to the density and the kinetic 
temperature of the gas. Therefore, with its observational advan- 
tage of a large number of lines being observable simultaneously, 
methanol is an excellent tool with which to probe the phy sical 
structure in dense molecular gas (also see Kama et al mm- 

A full spectral scan toward Orion KL covering fre- 
quency ranges 479.8 to 560.0 GHz (band la) and 1047.0 
tol 121.5 G Hz (band 4b) has been carried out by the HIFI 
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Guaranteed Time Key Program Herschel/HIFI Observations of 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



Extraordinary Sources: The Orion and Sagittarius B2 Star- 
forming Regions (HEXOS). Orion KL, located within the Orion 
molecu lar cloud at ~ 450 pc, is t he nearest massive star-forming 
region taenzel & Stutzki 1989h. It contai ns sever al kinematic 
components (e.g.. iBlake et aL Tl987; Perss on et al.l l2007). a hot 
core, multiple outflows, and the compact ridge, surrounding by 
cold quiescent gas. We have detected two methanol bands (i.e., 
groups of methanol lines within small frequency ranges) with 
transitions spanning a large range in upper state energy (80-900 
K) falling within a single HIFI bandpass (4GHz). In this paper, 
we show that these two bands of methanol can be used to de- 
termine the gas temperature and we explore thermal gradients in 
the Orion compact ridge. In §2, we discuss the HEXOS observa- 
tions and the analysis tools used in this study. In §3, we present 
the spectra of two methanol bands and their '^CHsOH counter- 
parts, and resolve the methanol emission from the compact ridge 
and the outflow. In §4, we show the population diagrams for the 
detected isolated methanol transitions for compact ridge and the 
outflow. And in §5, we examine how these two methanol bands 
can be used to provide additional information on the physical 
structure of the emitting regions. 



2. Observations 

The HIFI observations presented here, which are part of HEXOS 
(Bergin et al. 2010), are full spectral scans in bands la and 4b, 
covering a frequency range of 479.8 to 560.0 GHz and 1047.0 to 
1 121.5 GHz, respectively, pointing towards the Orion Hot Core 
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Fig. 1. HIFI spectra of two methanol bands in 
band la (521.7-525.6 GHz; upper panel) and 4b 
(1056.1-1065.6 GHz; lower panel). Red num- 
ber 1-42 label isolated methanol transitions, 
which are also listed in Table 1. Blue texts la- 
bel transitions from other molecules and the 
methanol transitions which are blended. Green 
texts and "B" label the methanol lines which 
are blended with different parity states. "U" la- 
bels two unidentified lines. 



at ay2000 = 5*35"' 14.3* and (5y2000 = -5°22'36.7". These ob- 
servations were carried out on March 1 and 3, 2010 using the 
dual beam switch (DBS) mode with redundancy 6, with the Wide 
Band Spectrometer (WBS) providing a spectral resolution of 1 . 1 
MHz over a 4 GHz IF bandwidth. The HIFI beam sizes are ~ 
43" in band la and 20" in band 4b, and the beam efficiency 
in both bands is 0.71. All data presented here were processed 
with HIPE ( Ott et alJl20"Tol) version 2.4, using the standard HIFI 
deconvolution {doDeconvolution task in HIPE). From line-free 
regions in the spectrum, we obtain a rms of T^* =0.065 K for 
band la and Ta*=0.116 K for band 4b. All spectra here were 
analyzed using the XCLASS program (https://www.astro.uni- 
koeln.de/projects/schilke/XCLASS), which provides access to 
the CDMS and JPL in olecular catalogs jMiiller et alj 1200 iL 
l2005t iPickett et al.l[T998l) for identi fying lines and obtaining the 
molec ular parameters of methanol dXu & Lovasll 19971 : IXu et al.l 
l2008h . 



3. Results 

3.1. Detection of distinctive metlianol bands 

There are 223 and 141 methanol lines with v,=0 and 1 de- 
tected in HIFI bands la and 4b, respectively. No molecular and 
frequency information for transitions with v,>l is available at 
present from the CDMS and JPL catalogs. Among these detected 
methanol lines, 112 in band la and 73 in band 4b are isolated 
transitions. The rest are transitions which are blended with either 
other methanol transitions or lines of other molecules. These de- 
tected isolated methanol transitions span a wide range in upper 
state energy (0-1000 K). The focus of this paper is to investigate 
how we can utilize the advantage offered by methanol - having 
many transitions within small frequency ranges - to probe the 
physical conditions in the emitting region. For this purpose, we 
look for groups of isolated methanol transitions which are con- 
fined within small frequency range in the spectrum but which 
also cover a wide range in upper state energy with minimal line 
blending. This approach is useful as it shows an example of how 
physical conditions can be probed by observing only a limited 
spectral range. 



Among all the detected methanol lines, there are two very 
distinctive methanol bands in HIFI bands la and 4b- both are 
constructed with a series of methanol lines confined within a 
few GHz range and have a similar characteristic spectral shape. 
Fig.[T]shows the spectra of these two methanol bands. The first 
covers 521.7 to 525.6 GHz in HIFI band la (upper panel). The 4 
GHz spectrum includes the methanol band peaked at 524.8 GHz, 
which is defined as the "524 GHz methanol band" in this paper. 
The second covers from 1056.1 to 1065.6 GHz in HIFI band 
4b (lower panel). The 10 GHz spectrum includes the methanol 
band peaked at 1061.6 GHz, which is defined as the "1061 GHz 
methanol band" in this paper. 

In the 524 GHz methanol band, 31 methanol lines are de- 
tected: 27 isolated methanol transitions (labeled in red numbers 
1-27), 2 lines either blended with other methanol transitions or 
transitions from other molecules (labeled in blue), and 2 lines 
which are blended with methanol transitions having the same 
quantum numbers (J, K, v,) but different parity with either E, 
A-I-, or A- state (defined as the blended "B" lines and labeled as 
"B" in green). Methanol has E and A symmetric states, which are 
considered as different species as no radiative transitions are al- 
lowed between them. A-type methanol can be further separated 
as + and - parity. Fig.[T]shows that the 524 GHz methanol band 
is dominated by isolated transitions, which are mostly AJ=0, 
K=-4— >-3 Q branch with v,=0 representing the torsional ground 
state. On the other hand, in the 1061 GHz methanol band, 42 
methanol lines are detected: 15 isolated transitions (labeled in 
red numbers 28-42), 6 lines blended with others, and 21 blended 
"B" lines. Fig.[T]shows that, unlike the 524 GHz methanol band, 
the 1061 GHz methanol band is dominated by the blended "B" 
lines mostly with AJ=0, K=7^6, v,=0 transitions. 

Table 1 lists the molecular parameters of the 42 isolated tran- 
sitions, including 37 E-type (33 v,=0 and 4 v,= l) and 5 A-type 
(4 v,=0 and 1 v,= l) methanol transitions, their peak intensities 
in units of antenna temperature (T^*), line frequency v, upper 
state energy E„, and S/u^ where S and fi are the line strength and 
dipole moment, respectively. The peak intensities are estimated 
with baseline continuum subtracted; the continuum level ~ 1 .23 
K for band la and ~ 6.45 K for band. These methanol lines have 
peak intensities between 0.4 and 12 K. They are widely spread 
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Fig. 2. Spectra of six isolated methanol lines. Each spectrum is 
plotted with the data (black) and the two-component Gaussian fit 
(red) to the data, which includes a narrower component (green) 
and wider component (blue). The number in the upper left corner 
and the energy (in K) in the upper right corner of each panel 
indicate the associated methanol transition and its upper state 
energy given in Table 1 . 



in upper state energy, which ranges from 80 to 900 K, covering 
almost the same energy range as those detected from the en- 
tirety of methanol lines detected in bands la and 4b. This energy 
range is sensitive to kinetic temperatures anticipated in molec- 
ular clouds, and these lines are thus a good set of transitions to 
probe the cloud conditions. We only focus on those 33 isolated 
E-type v,=0 (25 in band la and 8 in band 4b), in order not to 
get too involved in complicated excitation issues, any potential 
A/E abundance variations, and any v,>0 transitions that can be 
excited through radiative pumping of the strong Nvi2 torsional 
bands. 

Furthermore, among all 112 isolated methanol lines detected 
in the whole of band la, there are in total 55 isolated E-type 
Vf=0 transitions. This means that 25 (almost 50%) out of the to- 
tal 55 lines in band la (covering 80 GHz) are located within the 
4 GHz bandwidth of this 524 methanol band. On the other hand, 
among the 73 isolated methanol lines detected in whole band 4b, 
there are in total 36 isolated E-type v,=0 transitions- only ~ 20% 
of them are located in the 1061 GHz methanol band. The fact 
that the 524 GHz methanol band is dominated by isolated transi- 
tions while the 1061 GHz methanol band is mostly constructed 
with blended "B" lines, suggests that the 524 GHz methanol 
band is more suitable for exploration of the physical conditions. 
Hereafter, when we use the terms "524 GHz methanol band" and 
"1061 GHz methanol band" we are referring only to those iso- 
lated E-type v,=0 transitions. 

3.2. Resolving different kinematic components in Orion KL 

Fig. |2] plots the spectra of six isolated methanol lines with dif- 
ferent upper state energies. All methanol lines show asymmetry 
in the line profiles, indicating the existence of multiple contri- 
butions from different kinematic components in Orion KL, as 
suggested previously. For example, a high re solution VLA. spec- 
trum of one methanol line toward Orion KL dWilson et al.lll989h 
shows one narrower (linewidth AV~2 km s~', and V/,sfi=7.6 km 
s"') and one wider (AV~7 km s"', and V is r~(3. 2-9. 2 km s"') 
component. Consistent with previous results, the detected iso- 
lated methanol lines in Table 1 can all be fit with two Gaussians 
(see Fig. [2]): (1) a narrow component at V/,ss=8.6 km s"' with 
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Fig. 3. Spectra of four isolated methanol lines (black) plotted 
with their '^CHsOH counterparts (red), which are scaled by a 
factor of two in order to be compared to methanol lines more 
clearly. The number in the upper left corner and the energy (in 
K) in the upper right corner of each panel indicate the associated 
methanol transition and its upper state energy given in Table 1 . 

Table 1. Detected isolated CH3OH Lines 



label 


CH3OH 


y 


E„ 








V"* A.V/ J 


(MHz) 








1 


9^ . ^ i9S -) n P 

-^-^-4,0 ^ — ^-3,0 ^ 


521728.127 


841.28 


22.047 


0.44 


2 


24 A n— j24 3 n F 


522046 360 


783.41 


21.209 


0.45 


3 


^4,U ^^^3,0 ^ 


522121.802 


200.44 


1.997 


2.58 


4 


23 .n— >23 tr. F 


522340 009 

Z. Z. .J V/ . V/ V/ V 


727.83 


20.360 


0.74 


5 


21 ^n— >21 in E 


522861 608 


623.60 


18.646 


1.05 




20 An^lO 5n F 


523092 856 


574.94 


17.773 


1.28 


7 


14 1 r,^13nn E 


523274 293 


248.94 


10.900 


4.70 


8 


19 ,„_»19 ,„ F 


523306 240 


528.60 


16.897 


1.44 


9 


24, n—>23-,n A+ 


523346 243 


703.52 


10.320 


0.51 


10 


18-4,0 ^18-3 E 


523502 932 


484.56 


16.008 


1.59 


11 


17-4,0— > 17-3,0 E 


523683.978 


442.84 


15.112 


1.86 


12 


243,0^234,0 A- 


523744.752 


745.74 


7.788 


0.35 


13 


16-4,0^16-3,0 E 


523850.320 


403.43 


14.206 


2.18 


14 


15-4,0— » 15-3,0 E 


524002.815 


366.32 


13.290 


2.64 


15 


14-4,0^14-3,0 E 


524142.239 


331.53 


12.360 


2.92 


16 


13-4,0— >13-3,o E 


524269.306 


299.07 


11.423 


3.64 


17 


12-4,0— » 12-3,0 E 


524384.667 


268.91 


10.468 


4.20 


18 


11-4.0^11-3.0 E 


524488.918 


241.07 


9.499 


4.67 


19 


10-4.0^10-3,0 E 


524582.604 


215.55 


8.511 


5.32 


20 


24-2,0^24-1.0 E 


524614.631 


724.08 


20.883 


0.62 


21 


9-4.0^9-3,0 E 


524666.219 


192.34 


7.501 


5.90 


22 


8-4.0— *8-3,o E 


524740.210 


171.46 


6.461 


6.44 


23 


7-4,0—^7-3,0 E 


524804.977 


152.90 


5.380 


6.51 


24 


5-4,0^5-3,0 E 


524908.214 


122.73 


3.018 


6.26 


25 


4-4,0^4-3.0 E 


524947.258 


111.12 


1.647 


5.18 


26 


4-3,0^5-2,0 E 


525055.729 


85.92 


0.560 


3.21 


27 


13i,o—» 12-2,0 E 


525315.119 


232.30 


0.169 


0.57 


28 


11-4,0-^10-3,0 E 


1056354.886 


241.08 


8.659 


10.58 


29 


18o,i^l7i., E 


1056466.663 


696.17 


8.313 


2.32 


30 


22-1.0— >21-i o E 


1056627.148 


590.67 


17.746 


3.53 


31 


6-3.0—^5-2,0 E 


1057117.722 


111.46 


5.541 


11.21 


32 


22-2 1— *21-2 1 E 


1057831.005 


949.39 


17.720 


1.19 


33 


9-7J -^8-6.1 E 


1058132.994 


627.54 


12.064 


4.12 


34 


222,1 ->2l2,i A+ 


1058598.551 


883.67 


17.603 


1.37 


35 


22o,i— >21o.i E 


1058639.968 


885.60 


17.797 


2.23 


36 


13-2,0^12-1,0 E 


1059858.944 


237.30 


6.157 


7.47 


37 


22-7.0— >21-7.o E 


1061448.476 


846.67 


16.095 


1.37 


38 


22-5.0—^21-50 E 


1062221.782 


710.72 


16.893 


2.44 


39 


223,0^213.0 A+ 


1062261.800 


636.61 


17.458 


3.06 


40 


222,0— >2l2,o E 


1062312.231 


609.56 


17.571 


3.41 


41 


222,0^212,0 A+ 


1063916.882 


625.04 


17.739 


3.27 


42 


I81.0— >17o.o E 


1065427.650 


417.93 


7.792 


4.51 
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Fig. 4. Population diagrams of detected 33 E- 
type v,=0 methanol transitions for (a) the com- 
pact ridge and (b) the outflow. Transitions in 
HIFI bands la and 4b are plotted in black and 
red, respectively. Transitions of '''CH3OH are 
plotted in blue. Straight lines are fitted in order 
to obtain temperature estimates assuming LTE. 



AV=2.7 km s and (2) one wide component at V£sr=6-9 km Table 2. Detected '■'CHsOH counterparts 
s"' withAV=4-ll kms"'. Therefore, we fitted all methanollines 
with fixed W^sr at 8.6 km s"' and AV at 2.7 km s"' for the nar- 
row component, and freely fitting the Vlsr and AV for the wide 
component. It is well known that there are distinct spatial com- 
ponents (hot core, outflow, and compact ridge) along the line of 
sight that can be distinguished b y their kine matic sig nature in the 
molec ular emission profile (e.g.. lBlake et a l. 1987; P ersson et alJ 
Based on this information, we assign the V/,ss=8.6 km 
s~' component to the compact ridge, and the km s"^ 

component to the outflow. 



3.3. 



^^CH^OH counterparts for those detected methanol 
transitions 



In order to estimate the opacities of these '^CHaOH transitions, 
we look for their '^CHaOH counterparts throughout the whole 
band la. 13 of these 25 methanol transitions in HIFI band la 
with E„ < 370 K have '^CHsOH counterparts detected (Table 
2). We did not search '^CHsOH for those transitions in 4b be- 
cause no molecular information of '^CH30H transitions with 
frequency above 1 THz are available in either the CDMS and 
JPL catalogs. Fig.[3]plots four selected methanol lines and their 
'^CHsOH counterparts. Following the same strategy to separate 
different kinematic components as shown in §3.2, we then obtain 
the spectra for the compact ridge and the outflow component of 
each '^CHsOH fine. By comparing CH3OH and '^CHsOH lines, 
optical depth (r) for compact ridge and outflow components can 
be estimated using 



7^12 _ r„j2(i 

~3 ~ 



r«,i3(i-e-^'') 



(1) 



where we have assumed that CH3OH and '■'CHiOH have the 
same excitation temperature, CH3OH is optically thick, and 
Ti3=Ti2/f, where f is the 13C to 12C abundance ratio ~ 60 (e.g., 
ISheffer etal.ll2007h . Tn and To are the observed peak inten- 
sity for CH3OH and '^CHjOH lines. Table 2 gives the estimated 
CH3OH T12 for the compact ridge and the outflow, and we see 
that they are mostly < 10. We also estimate an upper limit of t\2 





"CHiOH 


V 


E„ 


Ta* 


T12 


Til 




(J^,v,) 


(MHz) 


(K) 


(K) 


(cr) 


(outflow) 


14 


-1,0— *13o,o E 


507380.068 


243.11 


0.58 


1.5 


12.9 


16- 


-4,0—^16-3,0 E 


528323.505 


396.25 


0.28 


16.9 


4.4 


15- 


-4,0^15-3,0 E 


528471.541 


360.05 


0.34 


7.8 


4.2 


14- 


-4,0— » 14-3,0 E 


528606.922 


326.10 


0.30 


7.2 


6.3 


13- 


-4,0—^13-3,0 E 


528730.329 


294.41 


0.41 


8.7 


4.8 


12- 


-4,0—^12-3,0 E 


528842.385 


264.98 


0.57 


6.1 


6.5 


11- 


-4,0^11-3,0 E 


528943.663 


237.82 


0.60 


5.6 


6.6 


9- 


-4,0^9-3,0 E 


529115.969 


190.27 


0.70 


5.7 


7.3 


8- 


-4,0— *8-3,o E 


529187.882 


169.89 


0.70 


4.5 


9.7 


5- 


-4,0— *5-3,o E 


529351.258 


122.33 


0.62 


5.3 


8.2 


4- 


-4,0— *4-3,o E 


529389.156 


111.00 


0.65 


6.4 


3.5 



~ 15 for those methanol lines with E„ > 370 K using the three 
times of the the band la rms (rms: T^i* =0.065 K) as the upper 
limit for all undetected '^CH30H emissions. 



4. Population diagrams of tlie isolated transitions 

Population diagrams are one traditional, but effective, way to dis- 
play the relative populations of different en ergy states, generally 
obtained from a set of observed transitions tooldsmith & Langej 
I1999h . This is a useful tool to study the physical conditions in 
the gas, as the relative populations of different states change as 
a function of kinetic temperature and density. A population dia- 
gram plots the upper state energy E„ as the x axis and the loga- 
rithm of N„/g„ as the y axis, where N„ and g„ are the column 
density and statistical weight in the upper state, respectively. 
N„ can be derived from the observed integrated line intensity of 
the associated transition with the eq. (9) in Goldsmith & Lange^ 
( 1999). Therefore, if all energy levels are thermalized with single 
kinetic temperature (LTE), the population diagram should be a 
straight line with the slope and y-intercept presenting the kinetic 
temperature and total column density. 

Fig.|4]shows the population diagrams, without correcting for 
possible beam dilution and opacity, for the detected isolated 33 
transitions of E-type v,=0 methanol and their '^CH30H counter- 
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Fig. 5. Population diagrams of the 524 GHz 
methanol band (in black) and the other 30 
isolated E-type methanol transitions detected 
throughout the rest of la (in red), for (a) the 
compact ridge and (b) the outflow, respectively. 
The fitted straight lines from Fig. |4] are also 
plotted for comparison. 



parts for the compact ridge (upper panel) and the outflow (lower 
panel), respectively. The straight lines shown in this figure are 
the slopes determined from least squares fitting with a linear 
function. We do not include the HIFI band 4b transitions in the 
fitting as they are observed with different beam sizes. This fig- 
ure shows that a curvature in the methanol population diagram is 
clearly seen for the compact ridge, which will be further exam- 
ined in §5. However, there are not enough '^CHsOH transitions 
to confirm whether this curvature is also shown for '^CHaOH in 
compact ridge. 

Unlike the curvature seen in the compact ridge data, the pop- 
ulations of the la transitions from the outflow can be satisfacto- 
rily represented by a single kinetic temperature in LTE of 1 84 K 
and 94 K for CH3OH and '^CHsOH, respectively. The '^CHjOH 
temperature might be a better approximation to the true kinetic 
temperature in this system because its emission is optically thin. 
However, the level populations derived from transitions within 
the 1061 GHz methanol band arising from the outflow exhibit a 
higher degree of scatter compared to those from compact ridge. 
This might suggest that the higher spatial resolution observations 
at 1 THz are less well coupled to the outflow and perhaps are 
contaminated by emission from other components such as the 
hot core. It is intriguing that the outflow in Orion KL is consis- 
tent with a single gas temperature, as outflow shocks are thought 
to be gas composed of a range of temperatures. It is likely that 
the derived temperature is simply showing an average gas tem- 
perature in this system. 

In addition, in order to verify whether the curvature in the 
compact ridge and the single temperature population in the out- 
flow are still valid for those transitions outside the selected dis- 
tinctive methanol bands, we compare the population diagram of 
the 524 GHz methanol band with that from the 30 isolated E- 
type v,=0 transitions in HIFI band la. Fig. |5] plots the popula- 
tion diagrams for the 524 GHz methanol band (in black) and 
the rest 30 lines (in red). A curvature is still present in com- 
pact ridge and all the level populations in the outflow can still 
be fitted with a straight line. The fitted straight lines indicating 
rotation temperatures of 125 and 184 K from Fig.|4]are plotted 
for comparison. This figure also shows that the addition of these 



30 transitions increases the scatter in the population diagram for 
both compact ridge and outflow. As these 30 transitions are ac- 
tually dominated by different Q branches, compared to those in 
the 524 GHz methanol band, this scatter is possibly due to a de- 
viation from LTE for populations within different K ladders. 

5. Probing the thermal structure in the Orion KL 
compact ridge 

As shown is §4, a curvature is present in the methanol popula- 
tion diagram of compact ridge, which cannot be fitted well with 
LTE excitation at a single temperature (Fig.|4]l. It has also been 
suggested in Fig. |5] that non-LTE excitation exists for methanol 
in the compact ridge, at least be tween transitions with different 
K ladders. Indeed. Leurini et al.l (|2004|) has shown that LTE can 
be a poor assumption for methanol even in clouds with high den- 
sities. However, to date, only transitions with up to upper state 
energies ~ 350 K have available collision rates for any non-LTE 
modeling; while the 524 GHz methanol band includes upper 
state energies from 80 to 900 K. In addition, according to Fig.|4] 
- 350 K is not a wide enough energy range to clearly distin- 
guish a curve from a straight line in the population diagram. 

Nevertheless, despite the currently limited capability for any 
non-LTE excitation modeling for all the detected transitions in 
the 524 GHz methanol band, it is intriguing that their popula- 
tion diagram shows such a cohesive curvature for compact ridge. 
This curvature might suggest presences of large opacity varia- 
tions, thermal gradients, and/or non-LTE excitation. Indeed, the 
opacity estimation in §3.3 (also see Table 2) is not negligible 
and can indeed affect the population diagram slope. In addition, 
thermal gradients are another potential contributor to this ob- 
served curvature. It has been suggested that the compact ridge is 
externally heated due to its very narrow linewidth and high tem- 
peratu re, inconsistent w ith normal centrally heated star-forming 
cores (lBlakeetalJll987h . 

In order to investigate how temperature gradients impact 
the curvature in a population diagram, we assume the com- 
pact ridge to be a spherical clump 5" - 10" in radius (r), hav- 
ing constant molecular hydrogen number density of 10^ cm"-*. 
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Fig. 6. Modeled spectra compared with the 
those observed in the 524 GHz methanol band 
(panel a), for models with constant tempera- 
ture of 50 (panel b), 100 (panel c), and 150 
K (panel d), assuming the compact ridge is 
a spherical clump 5" in radius (r), constant 
molecular hydrogen number density 10' cm"', 
and methanol abundance 10"*. 
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Fig. 7. Population diagrams of the observed 524 GHz methanol 
band (in black) and those models (in blue) shown in Fig.|6] 



and a constant methanol frac tional abundance of 10 ^ - 10 ^, 
as pre viously suggested (e.g. jBeuther et al.ll2006l : iPersson et all 
We consider various radial temperature profiles, and at 
each given point in the modeled clump, all energy levels are as- 
sumed to be populated in LTE with the temperature appropriate 
for that radius. With given level populations thus determined, 
we then use the ray propagation part of ID RATRAN program 
(iHogerheijde & van der Takl l2000) to generate a synthetic spec- 
trum for the 524 GHz methanol band observed with a 40" beam, 
consistent with the beam size in HIFI band la. 

Fig. |6] and Fig. |7] plot models with constant temperatures, 
to illustrate how the shape of this 524 GHz methanol band and 
its population diagram varies with different temperatures. Each 
modeled spectrum is constructed by combining the modeled 
compact ridge spectrum with the observed outflow spectrum. 
These two plots show that the observed spectrum cannot be ex- 
plained by a constant temperature in LTE. 

We next explore two approaches with nonuniform temper- 
atures: central heating and external heating. First, we examine 



whether a centrally-heated source can result in the observed 
spectrum with a temperature range from a few hundred K in the 
center of the clump decreasing to 30-50 K in its outer layer, 
which is cons istent with the quie scent gas temperature along the 
line of sight ( Bergin et alJll994l) . We allow any methanol abun- 
dance variation within the range 10"^ to 10"'. We were unable 
to find any models within these parameter ranges that provided 
a satisfactory match to the observed spectrum. The cold layer 
in the outer part of the clump always results in large opacities 
which produce much optically thicker lines than are observed. 

Second, we considered the possibility of externally-heated 
clumps and obtained a good fit that can reproduce the observed 
spectrum of the 524 GHz methanol band. The best fit is to 
a model of the compact ridge as a 7.5" (0.016 pc) spherical 
clump with a constant molecular hydrogen number density of 
10' cm"^, a temperature profile with constant 30 K for r<5.0", 
T~r^^ outside r>5.0", and a constant methanol abundance of 3 
X 10"^ and 3 x 10"' in the inner and outer part of the clump. 
This corresponds to a methanol column density of 9.5 x 10'^ 
cm"^ and 9.5 x 10" cm"^ in the inner and outer part of the 
clump, respectively. This gives a temperature of 120 K on the 
clump surface. We stress that the assumed temperature profile 
in this model is not based on any self-consistent calculation, 
but rather is a way to explore the range of gradients needed 
to obtain the observed level populations. PDR models do show 
that gas heating can prodti ce such sharp gradients as modelled 
here dKaufman et aljl 19991) and in the future we will examine 
more physically-motivated solutions. Fig.[8]plots the 524 GHz 
methanol band spectrum and the modeled spectrum for this best 
fit externally-heated solution in panel (a) and (b), respectively, 
and their associated population diagrams in Fig.|9^. We also ap- 
ply this model for transitions detected in HIFI band 4b to gener- 
ate its synthetic spectrum appropriate for a 20" beam, and fur- 
ther produce its population diagram shown in Fig.|9^. Fig. [8]: and 
Fig.|9}5 plot the modeled spectrum and its associated population 
diagram for one of the internally heated models we tested, for 
comparison, in order to show that the internally heated models 
can not explain the observed spectrum. 
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We conclude that the compact ridge is externally heated with 
the warm clump surface in front of the cold gas along the line 
of sight in LTE. Although our results show an external heated 
cloud can indeed result in a curvature in the observed population 
diagram, non-LTE excitation, such as far-infrared pumping, is 
also a possible contributor, which could potentially improve the 
fit for transitions having higher upper state energies. 

6. Conclusions 

We have examined the methanol emission in HIFI bands la and 
4b from our spectral scan observations toward Orion KL. In 
bands la and 4b, 223 and 141 methanol lines are identified. In 
band la, 112 out of the total 223 hues, and in band 4b, 73 out 
of the 141 Unes are isolated transitions. In order to minimize 
the confusion due to any blends, abundance variation between 
E- and A-type methanol, and optically pumping effects for v,>0 
transitions, we have focused on the isolated E-type v,=0 tran- 
sitions within two methanol bands: Methanol is known to have 
numerous transitions from far-infrared to millimeter regime, of- 
ten including different symmetric states (A and E), parities, and 
torsional slates, and blends. This makes the modeling and in- 
terpretation of underlying physical conditions difficult. But with 
the approach in this paper - focusing on one type of isolated 
methanol transitions - methanol can be a very useful tool to 
probe the physical conditions in molecular clouds. In particu- 
lar, the 524 GHz methanol band, including 25 isolated E-type 
V;=0 transitions with upper state energy ranging from 80 to 900 
K, provides a valuable set of lines that can be easily observed 
simultaneously and used to examine the temperature and other 
physical conditions. Furthermore, by resolving emission from 
different kinematic components (compact ridge and outflow) for 
each transition, we find a curved slope in the population diagram 
for the methanol emission from the compact ridge. We suggest 
that this can be explained by externally heating. Unlike the com- 
pact ridge, the population diagram of the outflow component 
shows a highly linear relation, suggesting a single temperature, 
but which might be the average temperature of the outflowing 
gas. Our study has shown that, with HIFI's broad instantaneous 
frequency coverage, these methanol bands can be very useful for 
proljing the structure of molecular clouds having various physi- 
cal conditions. 
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Fig. 8. (a) The 524 GHz methanol band spec- 
trum, (b) The best fit externally heated mod- 
eled spectrum, (c) Modeled spectrum of a cen- 
trally heated clump with temperature decreas- 
ing from 140 K in the center to 30 K on the 
clump surface. 
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Fig. 9. (a) Population diagrams for the HIFI 
band la (the 524 GHz methanol band, in black) 
and HIFI band 4b ( the 1061 GHz methanol 
band, in red) methanol bands and the best fit 
externally heated model for HIFI band la tran- 
sitions (in blue) and HIFI band 4b transitions 
(in green), (b) Population diagrams for the HIFI 
band la (the 524 GHz methanol band, in black) 
and HIFI band 4b ( the 1061 GHz methanol 
band, in red) methanol bands and the centrally 
heated model (la in blue and 4b in green) for a 
clump with temperature decreasing from 140 K 
in the center to 30 K on the clump surface. 



